Interference and Switching of Josephson current carried by nonlocal spin-entangled 
electrons in a SQUID-like system with quantum dots 
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Josephson current of spin-entangled electrons through the two branches of a SQUID-like structure 
with two quantum dots exhibits a magnetic-flux response different from the conventional Josephson 
current. Due to their interference, the period of maximum Josephson current changes from h/2e to 
h/e, which can be used for detecting the Cooper-pair splitting efficiency. The nonlocal spin entan- 
glement provides a quantum mechanical functionale for switching on and off this novel Josephson 
current, and explicitly a switch is formulated by including a pilot junction. It is shown that the 
device can be used to measure the magnitude of split-tunneling Josephson current. 
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Introduction - Nonlocal quantum entanglement becomes 
the focus of many investigations in recent years. The so- 
called Einstein-Podolsky-Rosen (EPR) pair— is the hall- 
mark of this phenomenon, which not only serves as the 
test case of violation of the Bell inequality 2 , but also 
works as the medium for quantum communication^ and 
quantum computation 4 . Photons are most intensively 
investigated, and their nonlocal entanglement has been 
successfully demonstrated^. It would be interesting to 
observe this phenomenon in electron systems, where en- 
tanglement may arise in either spin or spatial degrees of 
freedom. However, the generation and detection of non- 
local entangled electrons in solid-state systems is still a 
challenge, since electrons interact with the macroscopic 
Fermi sea around them and it is hard to control a par- 
ticular pair. 

A possible approach to achieve nonlocal entanglement 
of an electron pair is to take advantage of the intrin- 
sically spin entangled Cooper pairs in superconductiv- 
ity. Based on this idea, three-terminal devices consist- 
ing of an s-wave superconductor coupled to two leads 
made of quantum dots (QDs), Luttinger liquid, and nor- 
mal Fermi liquid have been theoretically proposed^— . 
It was shown that the Coulomb blockade effect may 
split a Cooper pair and force spin-entangled electrons 
tunnel into different leads. In a SQUID-like structure 
these spin-entangled electrons generate a new contribu- 
tion to Josephson current^. The anticipated Cooper pair 
splitting has then been explored experimentally in meso- 
scopic systems 1 - 1 - - — . The correlation between the resis- 
tances in the two leads due to the crossed Andreev effect 
(CAR) was used to measure the efficiency of Cooper pair 
splittingl^^ 2 .. 

Generally speaking, however, the resistance correlation 
does not provide a direct evidence for the nonlocal spin 
entanglement 7 -, since processes without spin entangle- 
ment may be involved 2 ^. In the present work, we pro- 
pose to detect nonlocal spin entanglement based on the 
interference of Josephson current, which is immune to 
any process involving disentangled electrons. For this 
purpose, we adopt the system in Fig. 1 which was first 
discussed by Ref. 10. We notice that the novel Joseph- 




FIG. 1: Schematic setup of a Cooper pair splitter. Two su- 
perconductors are connected by two leads with quantum dots 
(QDs) embedded. The QDs are defined by voltage gate gi 
and g2. 



son current carried by spin-entangled electrons tunneling 
through the two paths in Fig. 1 responds to the magnetic 
flux differently from the conventional Josephson current 
for which a Cooper pair tunnels through one of the two 
leads. Due to their interference the maximum Josephson 
current exhibits a variation of period h/e responding to 
the magnetic flux, in contrast of the conventional h/2e 
known for SQUID, which can be used to detect the split- 
ting efficiency of Cooper pair. The nonlocal spin entan- 
glement is then shown to provide a quantum mechanical 
functionale for switching on and off this split tunneling 
process. The function of a switch based on a pilot Joseph- 
son junction is formulated explicitly. It is shown that the 
device can be used to measure directly the magnitude of 
the novel Josephson current. 

Interference in presence of magnetic flux - It has been 
demonstrated experimentally-^ that QD can work as a 
notch to control the Cooper pair splitting, since the 
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Coulomb blockade effect suppresses the tunneling of two 
electrons of a Cooper pair through the same QD. The 
system can be described with a tunneling Hamiltonian 23 
H = H Q + Ht, with H = H$ + -Hqd, where Hs = 
Hl + Hr is the BCS Hamiltonian of the two supercon- 
ductors, 



with a = L,R for the left and right superconductor, a =t 
, 4 for the electron spin; Hqu = H u + Hd is the Anderson- 
type Hamiltonian for the QDs with one localized spin- 
degenerate energy level 



H V = e V X! a \o a r)<J +Un 



(2) 



and 5^* (uj) =<C a^^;a^'^ ^> is the retarded anoma- 
lous Green function of the QD. The superconductor is 
assumed to be a macroscopic system and thus its anoma- 
lous Green function is of the BCS form. The anomalous 
Matsubara Green function of the QD can be calculated 
with a contour integral method, 
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where Z is the partition function and /3 is the inverse tem- 
perature. Here we consider the regime e u = = e > 
where the ground state of the QD is empty of electron, 
thus avoiding the complexity induced by the local elec- 
trons on the QDs. Evaluating the Green function by per- 
turbation expansion in respect of the tunneling Hamilto- 
nian, we arrive at the lowest order result of the Josephson 
current at low temperatures 



with rj — u,d for the up and down QDs. In the presence 
of magnetic flux, the tunneling Hamiltonian is given by 
Ht = Htu + HTd, 

Htu = e^'^Ta^CR^fr^R) + e^^Tcl a (r u ^)a ua + h.c. 



7]-(j) 7]-$ 

h = [h sin(0 s - — ) + h sin(0 s + — ) + I 2 sin0 s ] (8) 
where <j) s is the superconducting phase difference, and 



H Td = e^o Ta da CR t(7 {r dtK ) + e^o Tcl a (r dih )a drT + h.c.,(3) 

with r^ jQ denoting the positions where QDs are con- 
nected to superconductors, T the tunneling matrix, $ 
the magnetic flux, and $o = h/2e the flux quantuam. 

The dc Josephson current can be evaluated using the 
Green function technique, and we arrive at the following 
three contributions Jj = I u + Id + lud- 
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with Ek = yf| + A 2 . It is diminished when the 
Coulomb interaction U and the superconducting energy 
gap A are large~, and 
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from Cooper pairs tunneling through the up and down 
QD respectively, which are similar to the result for a sin- 
gle QD Josephson junction 2 ^ except the Peierls factors, 
and 



lud = — —Re 



du; 

— n F {(j) 
Zn 
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fret 
ud 



( W )3 R et (^,Jr R )] 



from Cooper pairs which are split and tunnel coherently 
through different QDs. Here np(cj) is the Fermi distri- 
bution function, 5r a — r u a — r<j iQ is the terminal dis- 
tance, 9 R e *(aj,r) =<^ c R ^(r); c R ,j.(0) 3> is the retarded 
anomalous Green function of the right superconductors 
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FIG. 2: (color online). Magnetic response of maximum dc 
Josephson current. The dash-double dotted, dashed, dash- 
dotted, dotted, solid line is for the case of splitting efficiency 
7 = 0, 0.25, 0.5, 0.75, 1 respectively. Each curve was normal- 
ized and shifted for clarity. 
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which is controlled by the distances between the termi- 
nals of QD channels 1 -^. Here we note that the divergence 
of I2 on e is caused by the perturbation treatment, which 
can be eliminated by proper renormalizatio n 24 ' 25 . 

There are three contributions in the Josephson current 
in Eq. (8). The first two terms come from Cooper pairs 
co-tunneling from the up and down QD respectively. The 
interference between these two terms results in a period 
of $0 i n the magnetic response of maximum Josephson 
current, identical to that of a conventional SQUID. The 
third term is unique for the present system, in which 
the two electrons of a Cooper pair are split and tun- 
nel through different QDs, and the Peierls factors can- 
cel each other exactly. As the result, the interference 
pattern of maximum dc Josephson current exhibits a 
period of 2$o- It is clear that the behavior of the to- 
tal Josephson current is significantly modulated by the 
process of Cooper pair splitting. In order to see this 
influence more clearly, we define the splitting efficiency 
7 = I2/2I1. In the regime of strong Coulomb interaction 
and single particle resonate tunneling with e <C A <C U, 
one has 7 « (4 ^)e e~ 2l5r /^ sin 2 (fc F <5r)/(fc F Sr) 2 , where 
Sr = <5i-l = Str, where Sr = 5r^ = Sr-&, £ is the co- 
herence length and kp is the Fermi wave vector. The- 
oretically, the factor sin 2 (fcpi5r)/(fcF(5r) 2 sets the com- 
mon condition for observing the present quantity and the 
cross Andreev reflections 2 ^. The variation of maximum 
Josephson current with the splitting efficiency is illus- 
trated in Fig. 2. It is clear that with increasing splitting 
efficiency from zero to unity, the pattern evolves gradu- 
ally, and the period changes from $0 to 2$o- It is noticed 
that the splitting efficiency can be larger than unity, as- 
sociated with a similar curve to 7 = 1, except that the 
minima do not reach zero. 

Switching of the novel Josephson current - In order to 
contribute to the split-tunneling Josephson current, two 
single electrons tunneling through the two paths should 
be in a spin singlet state. This strong spin entanglement 
can be used for switching on and off this novel Josephson 
current. Let us set — [7/2 << —A < e < such that 
there is a localized electron on each QD 10 . For the split- 
tunneling Josephson current, it is the "on" /"off" state 
when the spins of the two localized electrons arc in a 
spin singlet /triplets state. The situation can be described 

byia 



/ = J 2 (i - S u • S d )sin< 



(11) 



Here we propose a way to control the two localized 
spins by introducing an additional pilot QD-junction as 
shown schematically in Fig. 3. The critical current J 
of the pilot junction varies continuously when the gate 
voltage on the QD is tune d 26 i 27 . The total Josephson 
current is then given by 



[2/ 1 +7 2 (i-S u -S d ) + J]sin<£ s (12) 




FIG. 3: (color online). Schematic setup of a switch of the 
novel Josephson current. The Cooper pair splitter displayed 
in Fig. 1 is connected in parallel with a pilot Josephson junc- 
tion. 



with the co-tunneling current 
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and the split-tunneling curren 
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where only the zero order terms in U are included for 
simplicity, since U is the largest energy scale. Both Ji 
and I2 are similar to the previous section, except that 
I\ is negative, which represents the 7r junction nature 
of each QD for co-tunneling current. For most cases, 
-h > h- 

The critical current of the total system can be straight- 
forwardly evaluated for the two spin configurations, 
which is plotted as function of the critical current J of 
the pilot junction in Fig. 4. For J < — 2Ii — I2/2, the 
critical current of the spin triplets state is larger than 
that of the singlet state, and vice versa. 

When a current is injected into the system through 
one of the superconductors, the superconductivity phase 
difference will be adjusted in order to pass this current 
without dissipation. If the current is small, the localized 
spins can take either singlet or triplet state. For cur- 
rent in between the two critical currents, the system will 
adjust the localized spins to achieve the larger critical 
current in order to reduce dissipation 10 . It is easy to fig- 
ure out that for J < —2I\ — I%/2 the spin triplets, thus 
"off", state is realized, and for J > — 2Jj — 1 2 /2 the spin 
singlet, thus "on" state, is realized, as depicted in Fig. 4. 

It is noticed that there is a gap I g in the maximal 
Josephson current of the total system when the critical 
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-2IfI 2 -2\ J 

FIG. 4: (color online). Maximum Josephson current of the to- 
tal system in response to critical current of the pilot junction 
J. The dashed/dash-dotted line is for the spin triplets/singlet 
state, and the solid lines are for the states realized. 



current of the pilot junction is swept, due to the split- 
tunneling Josephson current. By measuring the gap cur- 
rent we can evaluate the magnitude of the split-tunneling 
Josephson current since li = 2I g . At the switching point, 
the Josephson current carried by co-tunneling processes 
is suppressed significantly down to and thus the 

splitting efficiency of Cooper pair is enhanced largely. 
Before ending this section, we notice that the pilot QD 



can be replaced by a conventional SQUID, for which the 
critical current can be controlled by a magnetic field. 

Summary - In summary, we analyze the tunneling pro- 
cesses in a SQUID-like device with a quantum dot embed- 
ded in each junction, where electrons with nonlocal spin 
entanglement are generated and tunnel through different 
junctions, which contribute to the Josephson current. In 
presence of a magnetic flux, they carry a zero Peierls 
phase due to the cancelation between the two tunnel- 
ing paths, in contrast to the tunneling of Cooper pairs 
through one of the two junction. As the splitting effi- 
ciency increases, the period of the magnetic-flux response 
of maximum Josephson current changes from $o = h/2e 
to 2$ - The nonlocal spin entanglement is then used to 
switch the split tunneling process. The "on" ("off') state 
of the switch corresponds to the spin singlet (triplet) of 
the two localized electrons in the quantum dots of the 
Cooper-pair splitter, which can be controlled by a pi- 
lot junction. It is shown that the device can be used to 
measure directly the magnitude of the Josephson current 
carried by single electrons with nonlocal spin entangle- 
ment. 
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